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Dynamo action in a fully helical Beltrami (ABC) flow is studied using both direct numerical 
simulations and subgrid modeling. Sufficient scale separation is given in order to allow for large- 
scale magnetic energy build-up. Growth of magnetic energy obtains down to a magnetic Prandtl 
number Pm = Rm/Rv close to 0.005, where Rv and Rm are the kinetic and magnetic Reynolds 
numbers. The critical magnetic Reynolds number for dynamo action R C M seems to saturate at values 
close to 20. Detailed studies of the dependence of the amplitude of the saturated magnetic energy 
with Pm are presented. In order to decrease Pm, numerical experiments are conducted with either 
Rv or Rm kept constant. In the former case, the ratio of magnetic to kinetic energy saturates to 
a value slightly below unity as Pm decreases. Examination of energy spectra and structures in real 
space both reveal that quenching of the velocity by the large-scale magnetic field takes place, with 
an inverse cascade of magnetic helicity and a force-free field at large scale in the saturated regime. 

PACS numbers: 47.65.-d; 47.27.E-; 91.25.Cw; 95.30.Qd 



I. INTRODUCTION 

In recent years the increase in computing power, as well 
as the development of subgrid models for magnetohydro- 
dynamic (MHD) turbulence 0, 0, H H H has allowed the 
study of a numerically almost unexplored territory: the 
regime of low magnetic Prandtl number (Pm = Rm/Rv, 
where Ry and Rm are respectively the kinetic and mag- 
netic Reynolds numbers). This MHD regime is of partic- 
ular importance since several astrophysical Q and geo- 
physical @, § problems are characterized by Pm < 1 , as 
for example in the liquid core of planets such as Earth, 
or in the convection zone of solar-type stars. Also, liquid 
metals (e.g., mercury, sodium, or gallium) used in the 
laboratory in attempting t o g enerate dynamo magnetic 
fields arc in this regime |E EE EL E2] • 

In recent publications [H, EH EE EE E3 , driven turbu- 
lent MHD dynamos were studied numerically, within the 
framework of rectangular periodic boundary conditions. 
As Pm is lowered at fixed viscosity the magnetofluid be- 
comes more resistive than it is viscous, and it was found 
that magnetic fields were harder to excite by the dynamo 
process because of the increased turbulence in the fluid. 
The principal result was in obtaining the dependence 
of the critical magnetic Reynolds number R C M with the 
magnetic Prandtl number. These studies were done for 
several settings, ranging from coherent helical [lTj and 
non-helical [lj, EH forcing, as well as for random forc- 
ing pE EH- While for coherent forcing an asymptotic 
regime was found at small values of Pm, the behavior of 
R° M for random forcing is still unclear (see also [IE EE] 
for theoretical arguments based on the Kazantsev model 
0). 

For coherent forcing such as the Taylor-Green vortex 
(that corresponds to several laboratory experiments us- 
ing two counter- rotating disks), the value of R C M was ob- 
served to increase by a factor larger than six before the 
asymptotic regime for small values of Pm was reached 
[14(. Although the precise value of R C M in experiments is 



expected to be modified by the presence of boundaries, 
it is of interest to study what properties of the forcing 
can modify and decrease its value. It is well known from 
theory [2l[ , two-point closure models [22 1, and direct nu- 
merical simulations (DNS) at P M = 10, SH that 
the presence of net helicity in the flow helps the dynamo 
and decreases the value of R C M . 

In [l7j dynamos with a helical forcing function were 
studied using the Roberts flow, but mechanical energy 
was injected at a wavenumbcr |k| = y/2, which left little 
room in the spectrum for any back-transfer of magnetic 
helicity as expected in the helical case [IEHEUt]]; indeed, 
|k| = 1 is the only possibility since the computations are 
done in a box of length 27r corresponding to a k = 1 
gravest mode. In this work, in contrast, we study the ef- 
fect of a fully helical Arn'old-Childress-Beltrami (ABC) 
forcing [28| with energy injected at a slightly smaller scale 
(note that the ABC forcing is related to the Roberts forc- 
ing, since it can be defined as a superposition of three 
Roberts flows). As a result of the intermediate scale forc- 
ing, some a-effect or inverse cascade of magnetic helicity 
can a priori develop and a magnetic field at large scales 
can grow. 

ABC flows and helical dynamos were explored in many 
different contexts in the literaturefsee e.g. (29| for a 
study close to Pm = 1, and [2E [31| for studies in the 
context of fast dynamo action). The main aim of the 
present work is to study the impact of helical flows at 
intermediate scales in the development of magnetic fields 
through dynamo action at Pm < 1. In this context, it is 
worth noting that some simulations of ABC dynamos in 
the low magnetic Prandtl number regime were discussed 
in Refs. [IE [32|, although no systematic exploration of 
the space of parameters was attempted. Also, Ref. [33[ 
presented some preliminary results for the kinematic dy- 
namo regime with ABC forcing. In this work we will 
focus on the study of the generation of large scale mag- 
netic fields and of the non-linear saturation regime. A 
similar study was recently conducted in Ref. [34| using 
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mean field theory [26|, |27| and shell models. To the best 
of our knowledge, ours is the first attempt to systemati- 
cally study the saturation values of the fields for helical 
flows at Pm < 1 in numerical simulations. 



II. DEFINITIONS AND METHODOLOGY 

In a familiar set of dimcnsionless ( "Alfvenic" ) units the 
equations of magnetohydrodynamics are: 



dv 

~dt 
dB 

~dt ' 



v-Vv = -VP+jxB + !/V 2 v + f, (1) 



v VB = B • Vv + ?yV 2 B, 



(2) 



with V ■ v = V • B = 0. Here, v is the velocity field, 
regarded as incompressible, and B is the magnetic field, 
related to the electric current density j by j = V x B. 

V is the pressure, obtained by solving the Poisson equa- 
tion that results from taking the divergence of Eq. ([1]) 
and using the incomprcssibility condition V • v = 0. The 
viscosity v and magnetic diffusivity rj define mechanical 
Reynolds numbers and magnetic Reynolds numbers re- 
spectively as Ry = LU/v and Rm = LU/rj. Here U 
is a typical turbulent flow speed (the r.m.s. velocity in 

the following sections, U = (u 2 ) 1 ^ 2 , with the brackets 
denoting spatial averaging) , and L is a length scale asso- 
ciated with spatial variations of the large-scale flow (the 
integral length scale of the flow). We can also define a 
Taylor based Reynolds number R\ = AU/is, where A is 
the Taylor lcngthscale. defined below. 

Some global quantities will appear repeatedly in the 
next sections. These are the total energy (the sum of the 
kinetic Ey and magnetic Em energies) E = Ey + Em = 
\ J(u 2 + B 2 ) dV, the magnetic helicity H M = J A ■ B dV 
(where A is the vector potential, defined such as B = 

V x A), and the kinetic helicity Hy = J v • u> dV (where 
u) = V x v is the vorticity). While E and Hm are ideal 
[y = f] = 0) quadratic invariants of the MHD equations, 
Hk is not. In practice, kinetic helicity in helical dy- 
namos is injected into the flow by the mechanical forcing 
f (e.g., by rotation and stratification in geophysical and 
astrophysical flows [2l|V 

Equations ((T|) and ^ are solved numerically using a 
parallel pseudospectral code, as described in Refs. [la 
Il7| . We impose rectangular periodic boundary conditions 
throughout, using a three-dimensional box of edge 27r. 
The integral and Taylor scales are defined respectively as 



L = 27r^fc- 1 |v(k)| 2 /$>(k)| 5 



A = 2n ^|v(k)| 2 /^fc 2 |v(k)| 2 ^ 



1/2 



(3) 



(4) 



where v(k) is the amplitude of the mode with wave vector 
k (k = |k|) in the Fourier transform of v. 



The external forcing function f in Eq. {T]) injects both 
kinetic energy and kinetic helicity. For f we use the ABC 
flow 

fABC = fo{[Bcos(k F y) + Csin(k F z)]x+ 
+ [A sin(kpx) + C cos(kpz)] y + 
+ [A cos(k F x) + B sm(k F y)] z} , (5) 

with A = 0.9, B = 1, C = 1.1 323, and k F = 3. The 
ABC flow is an eigenfunction of the curl with eigenvalue 
k Fl and as a result if used as an initial condition it is 
an exact solution of the Euler equations. In the hydro- 
dynamic simulations, for large enough v (small Ry) the 
laminar solution is stable. As v is decreased the lami- 
nar flow becomes unstable and develops turbulence (see 
[35l | for a study of the early bifurcations at intermediate 
Reynolds numbers). 

To properly resolve the turbulent flow, the maximum 
wavenumber in the code k max = N/3 (N is the lin- 
ear resolution and the standard 2/3-rule for dealiasing 
is used) has to be smaller than the mechanic dissipation 
wavenumber k v = (e/z/ 3 ) 1 / 4 (e <~ U 3 /L is the energy 
injection rate). As a result, as v decreases and Ry in- 
creases, the linear resolution TV has to be increased. At 
some point the use of DNS to solve Eqs. (fTJ) and @ 
turns to be too expensive from the computational point 
of view and some kind of model for unresolved scales is 
needed. 

To extend the range of Ry and Pm studied, we use 
the Lagrangian average MHD e quat ions (LAMHD, also 
known as the MHD a-model) 0,111, 



cTv 
~dt 



Vv 



-VjVu{ — VP + j x B s 



-zA7 2 v 



dB s 

dt 



u s • VB S = B s • Vu s + ?yV 2 B. 



(G) 
(7) 



In these equations, the pressure V is determined, as be- 
fore, from the divergence of Eq. ([6]) and the incom- 
prcssibility condition. The subindex s denotes smoothed 
fields, related to the unsmoothed fields by 



v 
B 



= (l-Q'^V 2 )^ 

« 2 R V 2 )B S 



(1 



(8) 
(9) 



The total energy in this system is given by E = Ey + 
Em = \ J( v ■ u s + B ■ B s ) dV; it is one of the ideal 
quadratic invariants of the LAMHD equations. Equiv- 
alcntly, the magnetic helicity invariant is now given by 
Hm = J A s • B s dV, where the smooth vector potential 
is defined such as B s = V x A s . The expression for the 
kinetic helicity is the same in MHD and LAMHD. 

The LAMHD equations are a regularization of the 
MHD equations, and as a result they allow for simula- 
tions of turbulent flows at a given Reynolds number us- 
ing a lower resolution than in DNS. This subgrid model 
was validated against DNS of MHD flows in m 
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TABLE I: Parameters for the simulations: kinematic viscosity 
v, Taylor Reynolds number Rx, mechanical Reynolds number 
Rv, range of values of the magnetic Reynolds number Rm for 
a given flow, linear resolution N, and value of the mechanic 
filter length otv (with ay/as — [Rm / Rv] 3 ^ 4 in each case). 
For direct simulations (no subgrid model), ctv need not be 
defined. Runs in set 6a have the same values of parameters 
than in set 6, but a subgrid model at a lower resolution was 
used. Run 9 has the same parameters as set 9b with Rm — 41 
but was done using a DNS. The lowest Pm achieved for this 
set of runs is ~ 0.005. 



Set 


V 


Rx 


R\ ' 


Rm 


N 




1 


0.2 


11 


11 


9-16 


64 




2 


0.1 


21 


23 


10-19 


64 




3 


4 x 10~ 2 


55 


71 


17-71 


64 




4 


9 x 10~ 3 


161 


240 


18-54 


64 




5 


4 x 10~ 3 


250 


450 


15-450 


128 




6 


2 x 10~ 3 


360 


820 


10-41 


256 




6a 


2 x 10~ 3 


290 


840 


10-41 


64 


0.1 


7 


1 x 10" 3 


340 


1700 


14-42 


128 


0.0625 


8 


6.2 x 10~ 4 


680 


2500 


39 


512 




9 


5 x 10 -4 


500 


3400 


41 


512 




9b 


5 x 10~ 4 


500 


3400 


14-42 


256 


0.03125 


10 


2.5 x 10" 4 


1100 


6200 


77 


1024 





previous studies of dynamo action at low Pm, the ratio 
of the two filtering scales ay and cub was set using the 
ratio of the kinetic and magnetic dissipation scales, i.e. 
ay ja.B = P^ 4 [HI- The value of ay depends on the 
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linear resolution and was adjusted to l/cxv ~ h 

In the next section, we describe the computations and 
the results for both the kinematic dynamo regime [where 
j x B is negligible in Eq. |(Tj)], and for full MHD (where 
the Lorentz force modifies the flow). The first step is to 
establish what are the thresholds in Rm at which dynamo 
behavior sets in as Rv is raised and Pm is decreased 
(Section IlII Ap . The procedure to do this is the following 
(see e.g., Rcf. 01 )• Fhst a hydrodynamic simulation 
at a given value of i?y is done. Then, a small and ran- 
dom (non-helical) magnetic field is introduced, and sev- 
eral simulations are done changing only the value of Rm ■ 
At a given Rm, the magnetic energy can either decay 
or grow exponentially. In each simulation, the growth 
rate a is then defined as a = <1\xi(Em) / dt. The critical 
magnetic Reynolds number R C M for the onset of dynamo 
action corresponds to a = 0, and in practice is obtained 
from a linear interpolation between the two points with 
respectively positive and negative a closest to zero. The 
growth rate a is typically expressed in units of the recip- 
rocal of the large-scale eddy turnover time T = L/U. 

Once the values of R C M for different values of Pm < 1 
have been found, simulations for Rm > R%[ are con- 
ducted for longer times (Section llll Bp . In this case, mag- 
netic fields are initially amplified exponentially, and then 
saturate due to the back reaction of the magnetic field 
on the flow. In helical flows, this saturation is accom- 
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FIG. 1: Critical magnetic Reynolds R C M as a function of P M : 
DNS (solid line) and LAMHD (dotted line). Note the satu- 
ration for P M < 0.02. 
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FIG. 2: Growth rates as a function of Rm. Each line corre- 
sponds to several simulations at constant Rv (fixed v), and 
each point in the line indicates the exponential growth (or de- 
cay) rate at a fixed value of .Rm . The point where each curve 
crosses a = gives the threshold R C M for dynamo instability. 
Points from DNS are connected with solid lines, and labels 
are: set 1 (x), set 2 (□), set 3 (A), set 4 (o), set 5 (*), and 
set 6 (+). Points from LAMHD simulations are connected 
with dotted lines: set 6a (+), set 7 (x), and set 9b (□). Note 
the accumulation of lines near Rm ~ 20. 



panied by the growth of magnetic fields in the largest 
scale available in the box. In this regime, we will study 
the maximum value attained by the magnetic energy as 
a function of Pm (Section IIII C| . as well as the amount 
of magnetic energy at scales larger than the forcing scale 
(Section llll Dp . Finally, Section 4 is the conclusion. 



III. SIMULATIONS AND RESULTS 

In order to obtain a systematic study of dynamo action 
for ABC forcing and Pm < 1, a suite of several simula- 
tions was conducted. Table U shows the parameters used 
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Time 



FIG. 3: Time history of the magnetic energy for runs in set 
6 (dashed line) and in set 6a (solid line) with Rm ~ 41 and 
Pm = 5 x 10~ 2 . The inset shows the time evolution of the 
magnetic energy after the nonlinear saturation, in linear scale. 
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FIG. 4: Time history of the magnetic energy for runs in set 
5 (constant Rv ~ 450). The magnetic Reynolds number in 
each run is Rm ~ 22 (dashed line), Rm ~ 45 (dash-dotted 
line), Rm ~ 180 (dotted line), and Rm ~ 450 (solid line). 



in the simulations. Note that when a range is invoked in 
the values of Rm, it indicates several runs were done with 
the same value of Ry but changing the value of Rm to 
span the range (typically three to five runs). The set of 
runs 6 and 6a have the same parameters {y, T), and r.m.s 
velocity), but while set 6 comprises DNS at resolutions 
of 256 3 grid points, in set 6a the spatial resolution is 64 3 
and the LAMHD equations were used in order to further 
the testing of the model. Similar considerations apply to 
run 9 and set 9b. 



A. Threshold for dynamo action 

Figure [1] summarizes the results of the study of the de- 
pendence of the threshold R C M as Pm is decreased. For 
values of Rm above the curve, dynamo action takes place 
and initially small magnetic fields arc amplified. Below 
the curve, Ohmic dissipation is too large to sustain a 
dynamo. Noteworthy is the qualitative similarity of the 
curve between the ABC flow and previous results using 
different mechanical forcings [3, El [T3] - Namely, an 
increase in R° M is observed as turbulence develops, and 
then an asymptotic regime is found in which the value 
of R C M is independent of Pm- Note that LAMHD simu- 
lations were used to extend the study for values of Pm 
smaller than what can be studied using DNS. Simula- 
tions at the same value of Pm were carried with the two 
methods to validate the results from the subgrid model 
(sets 6 and 6a). This procedure was used before in Ref. 
jl4| . As in the previous study, the LAMHD equations 
slightly overestimate the value of R M . 

Besides the similarities in the shape of the curves for 
different forcing functions, two quantitative differences 
arc striking: (i) only a mild rise in R C M is observed here 
as Pm is decreased (a factor 2, while a factor larger than 
6 obtains for the Taylor-Green vortex (HI); and (ii) the 
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FIG. 5: Time history of (minus) the magnetic helicity for runs 
in set 5 (constant Rv ~ 450). The magnetic Reynolds number 
in each run is Rm ~ 22 (dashed line), Rm ~ 45 (dash-dotted 
line), Rm ~ 180 (dotted line), and Rm ~ 450 (solid line). 



asymptotic value of R C M for small values of Pm is ten 
times smaller than for other flows studied [ID, 53] ■ A 
similar result was obtained using mean field theory and 
shell models in Ref. [13], and the quantitative differ- 
ences observed were associated with the relative ease to 
excite large-scale helical dynamos compared with non- 
helical and small-scale dynamos. 

Note that the curve in Fig. [T]was constructed using the 
sets 1 to 7 and 9b of Table 1 . Several runs at constant Rv 
but varying Rm are required to define R C M . Set 9b reveals 
a dynamo at the lowest magnetic Prandtl number known 
today in numerical simulations, namely Pm = 4.7 x 10~ 3 . 

Figure [5] shows the details of how the thresholds for 
the determination of the R C M = f{P M l ) curve were cal- 
culated. For small initial Em, broadly distributed over 
a set of wavenumbers, r\ was decreased in steps to raise 
Rm in the same mechanical setting until a value of a w 
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was identified. A linear fit between the two points with 
a closest to provides a single point on the curves in 
Fig. [TJ Note that Figure [2] also gives bounds for the 
uncertainties in the determination of the threshold R C M 
(see e.g. Ref. errors in Fig. Q] can be defined as the 

distance between the value of R C M and the value of Rm in 
the simulation with a closest to 0. Note also the asymp- 
totic approach to a growth rate of order unity for large 
values of the magnetic Reynolds number, as for example 
in the runs in set 5. 



B. Time evolution 

A comparison of the time evolution of the magnetic 
energy in two dynamo runs with the same mechanic and 
magnetic Prandtl number (Rm ~ 41, Pm = 5 x 10~ 2 ) 
is shown in Fig. [3] One of the runs is a DNS from set 
6, while the other is a LAMHD simulation from set 6a. 
Two different stages can be identified at first sight in 
these runs. The kinematic regime at early times, with an 
exponential amplification of the magnetic energy (used 
to define the growth rates and thresholds in Figs. [T]and 
[2]), and the nonlinear saturated regime at late times. As 
expected from the results discussed in the previous sub- 
section, the LAMHD equations at a coarser grid (64 3 ) 
are able to capture the kinematic dynamo regime. While 
in the DNS with a resolution of 256 3 the growth rate is 
cr w 0.18, in the LAMHD simulation a w 0.20. But the 
LAMHD simulation also captures properly the nonlin- 
ear saturation (albeit the saturated level is reached a bit 
earlier) and the amplitudes of the magnetic energy in the 
steady state are comparable (see insert in Fig. |3|). Small 
differences observed in the time evolution are likely due 
to differences in the initial random magnetic seed. In 
the following, we shall use both DNS and LAMHD sim- 
ulations to study the nonlinear saturated regime at low 
Pm- 

In helical flows, as magnetic energy saturates, a large 
scale magnetic field develops (i.e., at scales larg er than 
the for cing scale) due to the helical a-effect |2ll . I22I [25], 
l2rl I27I . l3Qt | - It is of interest to know what happens with 
the amplitude of the magnetic field as the value of Pm 
is decreased. An example is shown in Figure 2J which 
gives the magnetic energy as a function of time for runs 
in set 5. Only the value of Rm (and therefore of Pm) 
is changed between the runs (Ry ~ 450 in all runs and 
Pm varying from 1 to 0.03). For large values of Rm 
(but not necessarily for values of Pm close to unity), the 
growth rate a is independent of Rm and of order one as 
noted in Section llll Al Furthermore, as Pm is decreased, 
both a and the saturation value of the magnetic energy 
decrease. However, for the lowest value of Pm, the mag- 
netic Reynolds number is quite low and in that context 
computations at a higher value of Ry and with the same 
sets of Pm are of value to see what fraction of the present 
result is a threshold effect at low Rm- 

Figure O shows the evolution of the magnetic helicity 
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FIG. 6: Time history of the magnetic energy in simulations 
at constant Rm ~ 41 (77 = 4 x 1CP 2 ). The different runs are 
taken from set 3 (Pm = 1, dash), set 4 (P M = 0.225, dot), 
set 5 (P M = 0.1, dash-dot), set 6 (P M = 0.05, solid), set 7 
(P u = 0.025, dash-triple dot), and finally set 9 (P M = 0.0125, 
long dash). 
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FIG. 7: Saturation value of the magnetic energy (normalized 
by the kinetic energy in the saturated regime) . The triangles 
correspond to simulations at constant Ry, while the squares 
correspond to simulations at constant Rm (squares connected 
with solid lines are from DNS, while squares connected with 
dotted lines are from LAMHD simulations) . Note the satura- 
tion at low Pm for constant Rm runs. 

as a function of time for the same simulations than in Fig. 
|4] The external forcing injects positive kinetic helicity in 
the flow. In the kinematic regime, the a effect is propor- 
tional to minus the kinetic helicity [26|, [2lJ ■ From mean 
field theory, the magnetic field in the large scales should 
grow with magnetic helicity of the same sign than the a 
effect (negative), as indeed observed (see Refs. (25l. |40]| 
for helical dynamo simulations at Pm = 1)- In the simu- 
lations, magnetic helicity grows exponentially during the 
kinematic regime. In runs with small Rm, the saturated 
state is reached shortly after the saturation of the expo- 
nential phase. But as Rm is increased, it is now clear 
that an intermediate stage develops in which magnetic 
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FIG. 8: Saturation value of the magnetic dissipation rate nor- 
malized by the kinetic energy in the saturated regime. The 
triangles correspond to simulations at constant Ry , while the 
squares correspond to simulations at constant Rm (squares 
connected with solid lines are from DNS, while squares con- 
nected with dotted lines are from LAMHD simulations). 



energy and helicity keep growing slowly. As a result, sat- 
uration takes place in longer times, and the time to reach 
the final steady state depends on the large scale magnetic 
diffusion time (T v « 47T 2 /?/). The dependence of the satu- 
ration time with Rm can be observed in Fig. [5] It is also 
worth mentioning that even in the runs with Pm < 1, 
the saturation of magnetic helicity can be well described 
by the formula ifjvf(i) ~ 1 — cxp[— 2r/fcQ(t — t sat )], where 
kg = 1 is the gravest mode, and t sat , is the saturation 
time of the small scale magnetic field [25| ■ This indicates 
that the slow saturation of the dynamo is dominated by 
the evolution of the magnetic helicity in the largest scale 
in the system. 

From Figs. Q]and[5]it seems apparent that small values 
of Pm have a negative impact on the amplitude of the 
magnetic field generated by the dynamo. However, dif- 
ferent results are obtained when the space of parameters 
is explored keeping Rm constant and increasing Ry , as 
another way to decrease Pm ■ Figure [6] shows the results 
in this case for the time evolution of the magnetic en- 
ergy. As Ry is increased from small values, a drop in the 
growth rate a and in the saturation value of the mag- 
netic energy is observed. But then an asymptotic regime 
is reached, in which both a and the saturation value seem 
to be roughly independent of Ry and Pm- As a result, 
we conclude that the behavior observed in Figs. 0] and [5] 
is the result of critical slowing down: if the space of pa- 
rameters is explored at constant i?y, as Pm is decreased 
Rm gets closer to R M until no dynamo action is possible. 
On the other hand, all the simulations with Pm < 0.05 
shown in Fig. [5] have Rm/R°m approximately constant 
(see Fig. [1} and critical slowing down is not observed. 
However, we note that the value of Rm for these runs is 
still modest. 
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FIG. 9: (Color online) Kinetic energy spectrum at t = [thick 
(blue) lines], and magnetic energy spectrum (thin lines) at 
different times: t = 11 (solid line), t = 29 (dashed line), 
t = 47 (dash-dotted line), t = 95 (dash-triple dotted line), 
t = 120 (long dashed line). The spectra are for a run in set 6 
with Rm ~ 41. The last time is in the saturation regime (see 
Fig. EJ. 



C. Saturation values 

The amplitude of the magnetic energy (normalized by 
the kinetic energy), after the nonlinear saturation takes 
place, as a function of the magnetic Prandtl number is 
shown in Fig. [7J This figure summarizes the results dis- 
cussed in Figs. [3]and[Sl As the value of Pm is decreased, 
if Ry is kept constant and Rm (and thus Pm) decreases, 
the saturation of the dynamo takes place for lower values 
of the magnetic energy. This is to be expected since as we 
decrease Pm we also decrease Rm and at some point R C M 
is reached. It is not clear whether such a strong depen- 
dence would be observed if the constant Ry runs were 
performed at substantially higher values of Ry as found 
in astronomical bodies and in the laboratory; however, 
such runs would be quite demanding from a numerical 
standpoint unless one resorts to LES (Large Eddy Sim- 
ulations) techniques, few of which have been developed 
in MHD (see e.g., 0, [H, H[). For values of Rm smaller 
than R M : 110 dynamo action is expected and the ratio 
Em/Ev should indeed go to zero. On the other hand, 
in the simulations with constant the ratio Em/Ev 
seems to saturate for P < 0.25 and reach an approxi- 
mately constant value close to sw 0.5. This indicates that 
small scale turbulent fluctuations in the velocity field are 
strongly quenched by the large scale magnetic field, as 
will be also shown later in the spectral evolution of the en- 
ergies. The ratio Em/Ev in helical large-scale dynamos 
is also expected to be dependent on the scale separation 
between the forcing wavenumber (here fixed to kp = 3) 
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FIG. 10: (Color online) Kinetic [thick (blue) line] and mag- 
netic energy spectra (thin line) at t — 210 in the simulation 
in set 6 with R M « 41 (P M = 0.05). The thick dashed line 
shows the total energy spectrum. 

and the largest wavenumber in the system (here k = 1). 
As the scale separation increases and there is more space 
for an inverse cascade of magnetic helicity, we expect the 
ratio Em/Ev in the Pm < 1 regime to also increase. 

Figure [5] shows the ratio of the magnetic energy dissi- 
pation rate cm = f] (j 2 ) to the kinetic energy dissipation 
rate ey = v (w 2 ) in the saturated state for the same runs 
than in Fig. [7] (in the LAMHD equations, the dissipa- 
tion rates are cm = f] (j 2 ) and ey = v (u> -u) s ), where 
w s = V x u s |). At constant Ry, for small values of 
Pm, critical slow down is again observed, as the value 
of Rm gets closer to the threshold. On the other hand, 
at constant Rm, more and more energy is dissipated by 
Ohmic dissipation as Pm is decreased. 



D. Spectral evolution 

In Refs. [l6|, it was shown using different forc- 
ing functions that even at low Pm the magnetic energy 
spectrum in the kinematic regime of the dynamo peaks 
at small scales. In these simulations, the critical mag- 
netic Reynolds number R C M was of the order of a few 
hundreds, and as a result small scales were excited. For 
ABC forcing, R C M is of the order of a few tens and close 
to the threshold small scales are damped fast. Only 
large-scale dynamo action is observed and thus, even at 
early times, the magnetic energy spectrum peaks at large 
scales. However, if Rm is increased above Rm ~ 400, a 
magnetic energy spectrum that peaks at scales smaller 
than the forcing scale (as in Refs. [l6|, [It}) is recovered. 
We focus here on large-scale dynamo action, and as a 
result will discuss the spectral evolution in simulations 




FIG. 11: (Color online) Kinetic [thick (blue) line] and mag- 
netic energy spectra (thin line) at t = 150 in the simulation in 
set 9 with Rm ~ 41 (Pm = 0.0125) in the saturated regime. 
The thick dashed line shows the total energy spectrum and 
the thin vertical dotted line the wavenumber at which the a 
filtering sets in. Note the compatibility of the spectra with a 
Kolmogorov law in the large scales for the kinetic spectrum, 
followed by a steeper power law. 

with Rm of a few tens. 

Figure [9] shows the evolution of the magnetic energy 
spectrum at different times for a run in set 6 with 
Rm ~ 41. As in previous studies, in the kinematic regime 
all the Fourier shells grow with the same rate. Then, 
magnetic saturation is reached and the mode at k = 1 
keeps growing until it eventually saturates itself. Figure 
[TOl shows the kinetic and magnetic energy spectra at late 
times (t — 210) after nonlinear saturation in the simula- 
tion in set 6 with R M tu 41 (P M = 0.05). At k = 1 the 
system is dominated by magnetic energy, but at smaller 
scales the magnetic energy spectrum drops fast. The ki- 
netic energy spectrum peaks at the forcing band (k = 3) 
and then drops with a slope compatible with k~ 3 . This 
drop is due to the action of the Lorentz force that re- 
moves mechanical energy from the k = 3 shell to sustain 
the magnetic field at k = 1 [U IHJ . 

A slope close to a fc -3 power law in the kinetic energy 
spectrum in the saturated regime at small scales is ob- 
served in several of the simulations with Pm < 1 ■ Simu- 
lations with small Pm and larger values of Ry were done 
using both the LAMHD equations and high-resolution 
DNS on grids of 512 3 and 1024 3 points (see Table H}. In 
these simulations, a power law close to fc~ 5 / 3 is observed 
before the kinetic energy spectrum drops to a steeper 
slope. As an example, Fig. QT] shows the kinetic and 
magnetic energy spectra in a simulation from set 9 us- 
ing the LAMHD model, with R M ~ 41 (P M = 0.0125). 
Slopes corresponding to fc -5 / 3 , k~ 3 , and fc~ 5 arc indi- 
cated as a reference in Fig. QT] A k~ 5 power law in 
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FIG. 12: Spectrum of relative magnetic helicity 
k~ Hmik)/ Ein(k) at different times in the simulation 
in set 6 with Rm ~ 41 (Pm = 0.05). The labels are as in 
Fig. [9] Note the evolution towards a force-free field at k = 1, 
the small excess of positive helicity at scales slightly smaller 
than the forcing scale, and the absence of relative magnetic 
helicity in the small scales at all times. 




FIG. 13: Spectrum of Hv{k) — k 2 HM(k), proportional to (mi- 
nus) the non-linear a-effect, in the simulation in set 6 with 
Rm » 41 (P M = 0.05). The labels are as in Figs. l9landlT2l 
the solid line is for t = 11, the dashed line for t — 29, and the 
long dashed line for t — 120. Note the drop of the spectrum 
at late times at scales smaller than fer- 
tile magnetic energy spectrum (following a fc -3 range) 
was observed in experiments of dynamo action with con- 
strained helical flows at low Rm [ill; i n addition, a fc~ 3 
power law for the kinetic energy spectrum is consistent 
with the observed magnetic energy spectrum [44]. Note 
that these power laws are only discussed here in order 
to be able to compare with the experimental data, but 
higher Reynolds numbers and thus more resolution will 
be needed in order to ascertain the spectral dependency 
of the flow and the magnetic field in the different inertial 
ranges of low Pm simulations. 

Figure [T2l shows the spectrum of relative magnetic he- 
licity k~ 1 HM(k)/EM{k) at different times for the same 



run as in Figs. [9] and [10] (run with Rm w 41 in set 
6). At all times, scales larger than the forcing scale have 
negative magnetic helicity, while scales of the order of, 
or smaller than the forcing scale have positive magnetic 
helicity. This is consistent with an inverse cascade of neg- 
ative magnetic helicity at wavenumbers smaller than /cf, 
and with a direct transfer of positive magnetic helicity 
at wavenumbers larger than fcp, as analyzed in [4oT | using 
transfer functions. The relative helicity in the k = 1 shell 
grows with time until reaching saturation. Note that at 
late times, HM(k = \)/EM(k = 1) ~ —1, indicating that 
the large scale magnetic field is nearly force-free. 

Figure Q2] also shows the spectrum of Hy(k) — 
k 2 HM(k), proportional to (minus) the non-linear a-effect 
[22| . Three times are shown for the same run than in 
Figs. El HOI and [TU (set 6, R M « 41). At early times 
(£ = 11 and t = 29) the spectrum of H v (k) - k 2 H M (k) 
is close to the spectrum of the kinetic helicity. However, 
as the large scale magnetic field grows (t = 120 is shown 
in the figure) the current helicity k 2 HM(k) quenches ki- 
netic helicity fluctuations and the total spectrum drops 
at scales smaller than kF- 

As a result, at late times the magnetic energy is mostly 
in the modes with wavenumber k = 1, which corresponds 
to the largest available scale in the system. In addition, 
the large scale magnetic field is force-free (maximum rel- 
ative helicity with HM(k) = EM(k) at k = 1. Figure [T4l 
shows slices of the velocity and magnetic fields at early 
and late times. The growth of a large scale magnetic 
field and the quenching of turbulent velocity fluctuations 
in the saturated regime can be easily identified. 

The situation resembles other inverse cascade situa- 
tions that have been studied numerically, in which the 
fundamental k = 1 mode dominates the dynamics at long 
times and its growth is only limited by its own dissipa- 
tion rate (2f| [3^, [4(| EB| • ln helical dynamo simulations 
at Pm = 1 this behavior has also been observed, although 
it was speculated that for Pm < 1 the inverse cascade of 
magnetic helicity and the generation of large scale fields 
should be quenched [25J- In fact, the generation of mag- 
netic energy at scales larger than the forcing scale is not 
quenched as Pm is decreased. This is illustrated in Fig. 
n~5T which shows the ratio of the magnetic energy in the 
k = 1 shell to the total kinetic energy in the saturated 
state as a function of Pm- Curves both at constant Rv 
and constant Rm are given. For constant Rm and small 
Pm the magnetic energy in the large scales seems to be 
independent of Pm and Rv- The overall shape of the 
curves is similar to the curves in Fig. indicating that 
at late times the evolution of the total magnetic energy 
is dominated by the magnetic field in the large scales. 



IV. CONCLUSION 

We have shown in this paper that the phenomenon 
of inverse cascade of magnetic helicity, and the ensuing 
growth of large-scale magnetic energy together with a 
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FIG. 14: (Color online) Plots of the velocity and magnetic fields in a cut at z — for the simulation in set 6 with Rm ~ 41 
and Ry ~ 820 (Pm = 0.05): (a) v z component in color and v x , v y indicated by arrows at early time, (b) same as in (a) for the 
magnetic field at early time, (c) same as in (a) at late time, and (d) same as in (b) at late time. 




1 10 100 

FIG. 15: Saturation value of the magnetic energy in the k — 1 
shell, normalized by the total kinetic energy. The triangles 
correspond to simulations at constant Rv, while the squares 
correspond to simulations at constant Rm (squares connected 
with solid lines are from DNS, while squares connected with 
dotted lines are from LAMHD simulations). 



force-free magnetic field at large times, is present at low 
magnetic Prandtl number, down to Pm — 0.005 in kine- 
matic regime studies and down to Pm = 0.01 in simu- 



lations up to the nonlinear saturation. The quenching 
of the velocity in the small scales, already observed in 
laboratory experiments, is also present. The augmenta- 
tion of the critical magnetic Reynolds number as Rv in- 
creases is less than in the non- helical case [III, [HI, HH , and 
even smaller than what was found for helical flows when 
the large-scale dynamo is not permitted, as e.g. for the 
Roberts flow at k « 1 [ljj • The reason for this difference 
is that in the present study we allowed for enough scale 
separation between the forcing scale and the largest scale 
for helical large-scale dynamo action to develop. The re- 
sults are in agreement with studies using mean-field the- 
ory and shell models to study both large- and small-scale 
dynamo action [34j . Large-scale helical dynamo action 
in the Pm < 1 regime requires much smaller magnetic 
Reynolds numbers to work than small-scale dynamos. 

The challenge remains, numerically, to be able to reach 
values of the magnetic Prandtl number comparable to 
those found in geophysics and astrophysics and in the 
laboratory, i.e. Pm ~ 10~ 5 . However, it is unlikely that 
the dynamo instability found here down to Pm = 0.005 
would disappear as Pm is lowered further. An open ques- 
tion, of importance from the experimental point of view 
when dealing with turbulent liquid metals, is whether the 
critical magnetic Reynolds number R C M will stabilize, for 
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a given flow, at a value intermediate between what it is at 
Pm = 1 and the peak of the curve (see Fig. [1]), or whether 
for large-scale helical dynamo action and extremely low 
values of Pm, it will go back down to the value it has at 
Pm = 1. The data up to this day suggests the former, 
but on the other hand the study made in the context of 
two-point closures of turbulence 44[ suggests the latter. 
This also means that reliable models of turbulent flows 
in MHD must be developed in order that we can explore 
in a more systematic way the parameter space character- 
istic of the flows of interest, as for the geo-dynamo or the 
solar dynamo. 
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